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We illustrate the use of effective theory techniques to describe processes involving unstable particles close to 
resonance. First, we present the main ideas in the context of a scalar resonance in an Abelian gauge- Yukawa 
model. We then outline the necessary modifications to describe VK-pair production close to threshold in electron- 
positron collisions. 



1. Introduction 

Processes involving heavy, unstable particles, 
play an important role in precision tests of the 
Standard Model and its hypothetical extensions. 
Typically these particles are studied in processes 
where they are produced close to resonance. It 
is thus important to overcome the difficulties re- 
lated to the breakdown of ordinary weak-coupling 
perturbation theory in the description of such 
processes. 

As is well known, the singularity of the inter- 
mediate propagator can be avoided if the width, 
r, of the unstable particle is taken into account 
through self-energy resummation. However, this 
alone does not guarantee an accurate description 
of the process. This is partly reflected by the 
fact that such a procedure may produce gauge- 
dependent results. Most approaches that have 
been put forward (for a summary see yQ) so 
far are focused on this particular aspect of the 
problem. Recently, we proposed a different ap- 
proach [2j making use of effective theory methods. 
This allows to perform a systematic expansion in 
the coupling as well as in T/M, where M is the 
mass of the unstable particle. The main advan- 
tage of this method is that it provides us with a 
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computational scheme for improving the accuracy 
of calculations in a systematic way, where gauge 
invariancc is automatic. 

2. Effective Theory 

At the heart of this method lies the observation 
that there is a hierarchy of scales between the 
mass and the width of the unstable particle. We 
exploit this hierarchy by constructing an effective 
theory and integrating out the modes that are not 
needed for the description of the external state. 
A first step to apply these methods in the context 
of unstable particles has been made in [3] and, as 
discussed in the next section, the programme has 
been carried out explicitly for a toy model in [2] . 

In technical terms, we identify the relevant 
modes and use the method of regions [2] to ex- 
pand the integrals in T/M. Together with the 
standard expansion in the coupling a we thus 
achieve a systematic organization of the calcu- 
lation in a series in a and T/M. 

As a first step in this procedure we inte- 
grate out hard momenta k ~ M. The effects 
of the hard momenta are incorporated into the 
matching coefficients, leaving the effective the- 
ory without dynamical hard modes. This is 
reminiscent and in fact often parallel to simi- 
lar developments in other effective theories such 
as non-relativistic QCD (NRQCD), heavy quark 
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effective theory (HQET) or soft-collinear effec- 
tive theory (SCET). The hard corrections are so 
cailed factorizable corrections, whereas the non- 
factorizable corrections are reproduced in the ef- 
fective theory by the still dynamical modes 
The precise nature of these modes depends to 
some extent on the underlying theory and even 
on the observable in question. In the next sec- 
tions we will give a more detailed discussion, first 
for a toy model, then for M^-pair production near 
threshold. 

3. An Abelian gauge- Yukawa model 

In order to outline the main steps in the con- 
struction of the effective theory (full details can 
be found in PJ) we consider a toy model, consist- 
ing of a heavy, scalar field <f> that decays through a 
Yukawa interaction into an "electron" and "neu- 
trino" . The electron field, tp, as well as the scalar 
field are charged under an Abelian gauge group. 
The Lagrangian is given by 



1 



M 2 
1 



A 



+ y ^x + y*^x^-^U) 2 + c ct , (1) 

where M and £ ct denote the renormalized mass 
and the counterterm Lagrangian respectively. 

For the purpose of illustration we outline the 
calculation of the totally inclusive cross section 
of the process D(q) e~(p) — > X as a function of 
s = (p + q) 2 with s - M 2 ~ ML ~ aM 2 , where 
a denotes collectively the gauge and Yukawa cou- 
plings. We obtain the cross section as the imag- 
inary part of the forward scattering amplitude 
T(s). The aim is to compute T through a sys- 
tematic expansion in a and 



M 2 



M 2 



r 

M 



(2) 



As a first step, we integrate out the hard modes 
with k ~ M. In the resulting effective theory, 
the effects of the hard modes are included in the 
matching coefficients of the operators. For the 
heavy scalar, this procedure is equivalent to the 



construction of HQET. In fact we write the mo- 
mentum of the scalar particle near resonance as 
= Mv>* + k^ where v 2 = 1 and k ~ M6 and 
define the soft field <p v (called "resonant" in 3 ) 
by removing the rapid spatial variation e r lMv%x 
from the </>-field. Thus, a soft field has momentum 
k ~ MS. We also define the matching coefficient 



A = 
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where s is the complex pole of the propagator. 
Since s and M are gauge invariant, A is guaran- 
teed to be gauge invariant as well. As required 
for a matching coefficient, A is given entirely by 
hard contributions. In fact, expanding the hard 
part of the self energy, Hh(s), as 

r2 gl n (*,0 



U h (s) = M 
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with I[( fc '') ~ a k , we can write the matching co 
efficient as 



(5) 
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The explicit form of the matching coefficient de- 
pends on the renormalization scheme. In the pole 
scheme A = — i T. 

The terms of the effective Lagrangian bilinear 
in (j) v can be written as 
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where we defined D^ T = — (v ■ D s )v fJ '. The 
interactions of the 0„-field with soft photons can 
be incorporated through the soft covariant deriva- 
tive D s = d — igA Sl because the separation into 
hard and soft parts respects gauge invariance. 

The Lagrangian Eq.® describes the propaga- 
tion of an unstable scalar particle close to res- 
onance. The first line is the leading term and 
the inclusion of A/2 into the propagator corre- 
sponds to self-energy resummation. The terms 
on the second line are suppressed by one power 
of S and contribute at NLO. Terms that are sup- 
pressed by further powers of S can be included 
systematically, if required. 
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In order to complete the construction of the 
effective Lagrangian we have to include kinetic 
terms for the soft and collinear photons and 
fermions, as well as the production and decay 
vertices for the unstable particle. In the present 
Yukawa model the NLO line shape is completely 
described by the addition of the production and 
decay vertices 

Ant = C [y<p v 1pn-Xn + + h.C.] 

)(Xn+i>n-), (7) 

where the matching coefficient C must be com- 
puted to one loop, and D = 1 at tree level. In 
this form of writing the effective Lagrangian we 
have integrated out the collinear modes, leaving 
only soft and external-collinear modes in the La- 
grangian. The latter describe soft fluctuations 
around the external light-like momenta. Finally, 
one must compute the scattering process to the 
one-loop order in the effective theory. We refer 
to for a thorough discussion. 

Outside the kinematic region 5 ~ a the ef- 
fective theory breaks down. To obtain a consis- 
tent description for all values of s, the result of 
the effective theory has to be matched to an off- 
resonance calculation in the full theory. 

4. PF-pair production close to threshold 

The toy model presented above is in many ways 
much simpler than the Standard Model, but it 
contains all the relevant features of resonant par- 
ticle production. Applying the same techniques 
to a realistic process may result in additional 
technical complications, but no new conceptional 
difficulties arise. 

An important application is W-pair production 
close to threshold at an electron-positron collider. 
This process is crucial for the precise determina- 
tion of the W mass and has been thoroughly stud- 
ied away from threshold. In particular, the one- 
loop electroweak corrections have been computed 
in the double pole approximation (DPA) .q] . Here 
we focus on obtaining results that are valid near 
threshold, where the DPA is supposed to break 
down. 

As mentioned before, the first step is to inte- 



grate out the hard modes. For pair production 
close to threshold, the W bosons are then de- 
scribed by a non-relativistic Lagrangian, analo- 
gous to NRQCD and the complications due to 
the Coulomb singularity can be addressed us- 
ing standard methods. We are thus left with 
dynamical degrees of freedom [1] familiar from 
NRQCD 0. We particularly emphasize the pres- 
ence of soft and potential gauge bosons. Working 
in the center of mass frame, the typical momen- 
tum of a W with velocity v is \k\ ~ Mw v and the 
typical (non-relativistic) energy is E ~ Mw v 2 . 
This corresponds to a potential mode. For soft 
modes, the energy and momentum scale as Mw v. 
The initial state fermions and the decay products 
are described by collinear modes familiar from 
SCET In fact we will need several collinear 
modes, one for each direction present in the pro- 
cess under consideration. 

To be specific, let us consider the process 

e-( Pl )e + (p 2 )^ ^-(h)^(l 2 )u(l 3 )d(U) (8) 

with {h+l 2 ?-M w ~ (lz+l A f-M w ~ M w v 2 ~ 
Mw^w- By assumption \fs/2 — Mw ~ Mw^w- 

At leading order, the description of the process 
can be split into three parts. First, the VL^-pair is 
produced, then it propagates, and finally the W 
bosons decay. At higher orders, this simple pic- 
ture is complicated by the exchange of Coulomb 
photons, of photons that connect the various 
stages and/or by the presence of additional pho- 
tons in the final state. Furthermore, single- 
resonant and non-resonant diagrams have to be 
taken into account. We count a s ~ v ~ y/a ew 
and by next-to-leading order (NLO) contributions 
we understand the 0(a ew ) ~ 0{v 2 ) ~ 0(a 2 ) cor- 
rections. In this article we present all contribu- 
tions that are needed at yNLO, that is, all cor- 
rections of order v ~ a s ■ 

Starting with the construction of the non- 
relativistic Lagrangian we first remark that the 
gauge dependence of the vector-boson propagator 
is not an issue. Indeed, expanding the momen- 
tum of a non-relativistic W as fc M = Mw f M + q^, 
where = (1,0), and q 11 = (q°,q) and taking 
into account q° ~ Mw v 2 and \q \ <~ Mw v the 
propagator in a general gauge upon expansion 
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The propagator scales as v~ 2 , is gauge inde- 
pendent and describes three polarization states 
of a non-relativistic particle. This is consis- 
tent with the fact that the degrees of freedom 
with mass \/£Mw have been integrated out. 
(In Feynman gauge the propagator is given by 
-ig»"/(2M w q° - q 2 ). The effects of the un- 
physical scalar polarization are cancelled by the 
pseudo-Goldstone field, which has mass My/ and 
is still present in the effective theory.) Including 
the decay width we obtain for the non-relativistic 
Lagrangian up to NLO 
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where Q± denote the non-relativistic vector fields 
with mass dimension 3/2 for the bosons. The 
matching coefficient A = (s — Ml) /My/ is de- 
fined as before. Its leading part, ~ Mwot ew 
scales as D° and D 2 and, thus, has to be included 
in the leading order Lagrangian. Accordingly, the 
propagator of a f2± with energy E and momen- 
tum k is given by 
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Higher order corrections to A can either be re- 
summed, i.e. included in the propagator, or in- 
cluded perturbatively as interactions. At VNLO 
we need A^ 3 / 2 ), the 0(a ew a s ) corrections to A. 
Higher order corrections are at least NLO as are 
the terms in the second line of Eq. (|10fl . 

The production of a W~ W + pair is described 
by effective vertices to be added to C nr . In order 
to obtain the corresponding operators and their 
matching coefficients, we have to compute the 
amputated, renormalized on-shell Green function 



for e + e~ — > W + W~ to the desired order in ordi- 
nary weak-coupling perturbation theory. At lead- 
ing order, only the helicity configuration e^e^ 
contributes and the corresponding operator reads 
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where (s-l) is an external-collinear field [2j with 
large momentum in the ft (—ft) direction and we 
introduced the notation a^b^ = a l V + a^b 1 . In- 
cluding terms that are suppressed by one power 
of v we get additional operators 

4 1/2) = $r i^ j e L ) (^(-i)Dm*j) (is) 
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with the matching coefficients 
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(14) 
(15) 

(16) 
(17) 

All derivatives on the potential fields in 

Eci. (|13H scale as D l ~ My/ v. We also remark that 
the e \ helicity configuration does not vanish 
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at VNLO. Consequently, there are additional op- 
erators in yCp 1 ^ involving en and e^. 

The decay of the vector bosons can be de- 
scribed in a similar way as the production, namely 
through decay vertices. The operators with 
their matching coefficients are obtained through 
matching of the corresponding on-shell Green 
function. For the leptonic decay, at VNLO there 
are no loop corrections to be taken into account, 
whereas for the hadronic decay, there are correc- 
tions of the order a s ~ v. When the calculation is 
completely inclusive on the hadronic decay prod- 
ucts, these corrections, together with the correc- 
tions from final state emission of an additional 
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gluon, are taken into account by the hadronic part 
of the decay width. For the decay operators that 
are relevant to the process Eq.© we obtain 

4 0) = + tf+urfdi) . (18) 

Kinematic corrections to C? are suppressed by 
at least two powers of v. 




Figure 1. Leading order Feynman diagram in the 
effective theory. 

We are now in a position to compute the lead- 
ing order amplitude within the effective theory. 
The corresponding diagram is depicted in Fig- 
ure^ For the vertices we have to take the leading 
order production and decay vertices, Eas. i|12|) and 
(|18|l . The propagators of the intermediate vector 
bosons are given by Ea. Hll|) . Putting everything 
together we obtain 

^ (0) = *c^_ ipa _ lnliryJl]pi _ )x (1Q) 

(h-\l l \h-)(h-W\h-) 

( rp (h+h) 2 aw\ ( rp (h+u) 2 awV 
\ ~2aT^ —) \ t *--5Kfo —) 

where E x = + 1% - M w , E 2 = l^ + ll~ M w and 
we used standard helicity notation. 

Turning to the calculation of the v'NLO ampli- 
tude, we also have to consider the helicity config- 
uration e^e\. However, the corresponding cal- 
culation is analogous to the calculation of A^°\ 
barring the replacement of the production ver- 
tex. Let us, therefore, focus on the corrections 
needed in the case of ej. 

For all stages, production, propagation and de- 
cay we have to include corrections. In the case 
of the production stage, we have to include dia- 
grams as shown in Figure ^ with the production 



vertex due to operators given in Ea. (|13|l rather 
than Ea.Q12[l. For the propagation stage, we need 
to include the 0(a ew a s ) corrections to A. If re- 
summed, this results in a change in the propaga- 
tor AW -> AW + A< 3 / 2 ). Finally, regarding the 
decay stage, we mention again that the 0(a s ) 
corrections to the decay vertex are not explicitly 
needed in a completely inclusive calculation. 

Apart from these trivial corrections, there is 
a further contribution at \/NLO, namely the ex- 
change of a single Coulomb (potential) photon, 
shown in Figure This contribution is sup- 
pressed by a/o ~ u relative to the leading order 
amplitude. Thus, contrary to top pair produc- 
tion close to threshold, these contributions can be 
treated perturbatively and need not be summed. 
This is the only correction at vNLO that is not 
due to hard modes. 




Figure 2. The exchange of a single Coulomb pho- 
ton is suppressed by a/v relative to the leading 
order amplitude and, thus, contributes at ^/NLO. 



The momentum of the potential Coulomb pho- 
ton scales as q° ~ A%» 2 , \q\ ~ Mwv. Thus 
the propagator is given by ijq 1 ■ Reading off the 
Feynman rules of the WW^ vertex of Ea. flOp . 
we obtain for the amplitude due to single photon 
exchange 

^(V2,c) _ _■ ( 4?ra ) A (0) x ( 20 ) 
r d D q 1 1 

1 

where k = kx = 1% + 1% and ki =^3 + ^4 = —k. 
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After performing the q° contour integral, the q 
integration is straightforward and we obtain 

A (l/2,c) = A (0) "M* x (21) 

1*1 

, \k\ 

arctan — - , 

y/M w (AW -E x - E 2 ) - it 

in agreement with [Sj. 

We note that it is possible to resum the con- 
tributions due to multiple potential photon ex- 
change and derive a Green function for the W 
pair. This is described most naturally within the 
context of a potential non-relativistic Lagrangian, 
analogous to potential non-relativistic QED 9 . 

5. Conclusions 

The fundamental reason for the breakdown of 
ordinary weak-coupling perturbation theory for 
processes involving resonant unstable particles is 
the appearance of a second small scale. We de- 
scribed a method that overcomes these problems 
using an effective theory approach and applied 
this first to a toy model involving an unstable 
scalar particle and then to VF-pair production 
near threshold. 

The applications we presented are technically 
rather simple. In the case of VF-pair production 
this is because we have considered only the 0(v) 
corrections. However, we would like to stress once 
more, that including higher-order corrections in- 
troduces only technical problems and no new con- 
ceptual problems arise. 

In this article, we have restricted ourselves to 
a totally inclusive cross section. More compli- 
cated final state kinematics requires in general 
the introduction of additional modes in the effec- 
tive theory. Also, phase-space integrals have to 
be expanded. This is done most conveniently by 
working directly with cut-diagrams, rather than 
amplitudes, because this facilitates the use of the 
same methods as for loop integrals. Thus, the 
method presented here is not restricted to a small 
number of special cases. It provides us with a 
consistent systematic computational scheme that 
can be applied to a wide variety of processes. 
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